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such as bangs, bounces, crunches, rips, sudden singularities, and extremality events. 
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the class of bangs, cruncjes, and sudden singularities for which the dominant energy 
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1. Introduction 

The cosmological principle, which is ultimately a distillation of our knowledge of 
observational cosmology, leads one to consider cosmological spacetimes of the idealized 
FRW form [1, 2, 3, 4]: 

ds 2 = -dt 2 + a{t) 2 l-^r-, + r 2 [d9 2 + sin 2 9 d0 2 ]l . (1) 
[ 1 — kr 2 J 

The central question in cosmology is now the prediction, or rather retro-diction, of 
the history of the scale factor a(i) [5, 6]. The key dynamical equations (assuming 
applicability of the Einstein equations of general relativity to cosmology in the large, 
an assumption which at least some cosmologists are beginning to question) are the 
Friedmann equations (which redundantly imply the conservation of stress-energy). 
Indeed, in units where SnG^ = 1 and c = 1, we have 

P(t) = 3 g + * ) , (2) 

P(* =-2---5--5, ( 3 ) 

a a 2 a 2 

p(t) + 3p(t) = -6-, (4) 
a 

and the related conservation equation 

p(t) a 3 + 3[p(t) +p(t)]a 2 a = 0, (5) 

where any two of equations (2, 3, 4) imply (5). § Until quite recently the only FRW 
cosmological singularities receiving significant attention were the big bang (and its time- 
reversed partner, the big crunch) where a(t) — > 0. Recently, the class of interesting 
cosmological singularities has become much broader [7, 8, 9, 10, 11, 12, 13, 14, 15]. 
In this article we shall pin down necessary and sufficient conditions for the various 
cosmological singularities and related events, (which we shall generically refer to as 
cosmological milestones), both at the level of kinematics, and at the level of dynamics. 

• The kinematical analysis boils down to developing precise definitions of these 
cosmological milestones ("major events") in terms of the scale factor a(t). 

• The dynamical analysis uses the Friedmann equations to analyze the classical 
energy conditions in the vicinity of these milestones in a model-independent manner 
(without even assuming the existence of any equation of state). 

§ If in contrast you try to use equations (5) and one of (2, 3, 4) to derive the other two equations 
you may need to make some specific choice of integration constants. For instance, (5) and (4) do not 
contain k, and in this situation k arises as an integration constant, which can without loss of generality 
be set to k = — 1/0/ + 1 by suitably rcscaling a(t). Secondly, and more significantly, combining (5) 
and (3) reproduces (2) only up to an ambiguity, p ~ K/a 3 , corresponding to an arbitrary quantity of 
"dust" which has to be eliminated "by hand" . 



Necessary and sufficient conditions for big bangs, bounces, crunches, rips. 



3 



The overarching aim of the article is to provide a (hopefully complete) catalogue of the 
interesting cosmological milestones that might be encountered in unusual cosmological 
models, and determine (with an absolute minimum of technical assumptions) the 
necessary and sufficient conditions for their occurrence. 

2. Kinematics — definitions 

Generic cosmological milestone: Suppose we have some unspecified generic cosmologi- 
cal milestone, that is defined in terms of the behaviour of the scale factor a(t), and which 
occurs at some finite time t . We will assume that in the vicinity of the milestone the 
scale factor has a (possibly one-sided) generalized power series expansion of the form 

a(t) = co\t - t r + d\t - t Q \ m + c 2 \t - t Q \ m + c 3 \t - t Q \ V3 + . . . (6) 

where the indicial exponents rji are generically real (but are often non-integer) and 
without loss of generality are ordered in such a way that they satisfy 

rj < m < rj 2 < rj 3 . . . (7) 

Finally we can also without loss of generality set 

c > 0. (8) 

There are no a priori constraints on the signs of the other q, though by definition q ^ 0. 

Physically, this generalized power series expansion is sufficient to encompass all the 
models we are aware of in the literature, and the indicial exponents rji will be used 
to classify the type of cosmological milestone we are dealing with. Indeed the lowest 
few of the indicial exponents are sufficient to determine the relationship between these 
cosmological milestones and the curvature singularities more traditionally encountered 
in the general relativity literature. Mathematically, power series of this type generalize 
the notions of Taylor series, Laurent series, and Frobenius series, and are even more 
general than the generalized Frobenius series adopted in [16]. The current definition is 
motivated in part by the fact that generalized Frobenius series (compositions of power- 
law prefactors with Taylor series in some variable raised to a possibly different exponent) 
commonly occur when expanding solutions of differential equations around their singular 
points. 

For many of our calculations the first term in the expansion is dominant, but even 
for the most subtle of the explicit calculations below it will be sufficient to keep only 
the first three terms of the expansion: 

a(t) = c \t - t^ + Cl \t - t Q \ m +c 2 \t-t Q \ m Vo<Vi<V2] c > 0. (9) 

Big bangs and big crunches: The most basic of the cosmological milestones are big 
bangs and big crunches, for which the scale factor a(t) — > at some finite time as we 
move to the past or future. Let the time of the big bang (if one occurs) be denoted by 
t* and the time of the big crunch (if one occurs) be denoted by t®. We shall say that 
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the bang or crunch behaves with indicial exponents (0 < 770 < Vi ■ ■ •) if the scale factor 
possesses a generalized power series in the vicinity of the singularity. That is, if 

a(t) = co(t - Q V0 + ci(f - U)^ + ... (10) 

or 

a(t) = c (t® - ty° + ci(t® - tr + . . . (11) 

respectively. Note that the series have been carefully constructed to make a(t*) = and 
a(i®) = 0. 

Big rips: A "big rip" is said to occur if a(t) — > 00 at finite time [7, 8]. We can 
distinguish a "future rip" from a "past rip", where the literature to date has solely 
considered future rips (as a past rip would be a most unusual and unexpected beginning 
to the history of the universe). Let the time of the rip, if it occurs, be denoted t#, then 
we define the indicial exponents of the rip (either future or past) to be (170 < 771 . . .) if 
the scale factor possesses a generalized power series in the vicinity of the rip: 

a(t) = c |t* - + ci|f« - t\ m + . . . , (12) 

with rjo < and Co 7^ 0. Note the similarity to bangs and crunches, with the only 
difference being in the sign of the exponent i] . Note that the series has been carefully 
constructed to make a(t&) = 00. 

Sudden singularities: Past or future sudden singularities are a recent addition to the 
collection of cosmological milestones, and are defined by some time derivative of the scale 
factor diverging at finite time, while the scale factor itself remains finite [11, 12, 13]. 
Again, almost all attention has been confined to future sudden singularities as (a past 
sudden singularity would be a most unusual and disturbing beginning to the history 
of the universe). Let the time of the sudden singularity, if one occurs, be t Q (past or 
future). A suitable definition of the exponent of a sudden singularity is to take r] = 
and rji > to give 

a(t) =c + Ci|f-f Q |" 1 + ... (13) 

with c > and 771 non-integer. Thus a(t Q ) = c is finite and a sufficient number of 
differentiations yields || 

a W(t _> t Q ) ~ c 7 7l (7 7l -l)( Vl -2)...( Vl -n + 1) \t - t r~ n -> 00. (15) 

II In particular the toy model considered by Barrow [11, 12, 13] can be written 

a(t) = c [(* e - ty - 1] + co(t - Uf (14) 

and falls into this classification when expanded around the time of the sudden singularity t Q , while it 
falls into the classification of big bang singularities considered above when expanded around the time 
of the big bang, t*. 
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Extremality events: Other common cosmological milestones are extremality events, 
which are not singularities in any sense. For these events a(t) exhibits a local extremum 
at finite time, so in particular a — + 0. In the vicinity of extremality events we can model 
the scale factor using ordinary Taylor series so that in terms of our generalized series 
we have i]q = and r]i G Z + . In particular 

• A "bounce" is any local minimum of a(t), the time of such an event being denoted 
by t„ so that a^(t.) = is zero [17, 18, 19, 20]. The "order" of the bounce is the 
first nonzero integer n for which the 2n'th time derivative is strictly positive [18, 19]: 

a^ 2n \t.) > 0, (16) 

so that 

ait) = ait.) + j^a^{t.) [t - t.] 2n + .... (17) 

• A "turnaround" is any local maximum of a(t), the time of such an event being 
denoted by t , so that a^(t Q ) = is zero [21]. The "order" of the turnaround is 
the first nonzero integer n for which the 2n'th time derivative is strictly negative: 

a (2n) (*o)<0, (18) 

so that 

a(t) = a(t ) + J-^ a ( 2n \t ) [t - t ] 2n + .... (19) 

• An "inflexion event" is an extremality event that is neither a local maximum or a 
local minimum, the time of such an event being denoted by t@. (Inflexion events 
can be thought of as an extreme case of "loitering" [22], in the limit where the 
Hubble parameter momentarily vanishes at the inflexion event.) The order of the 
inflexion event is the first nonzero n for which 

a (2n+1) (* @ ) ± 0, (20) 

so that 

a(t) = a(t @ ) + 7^-^ (2n+1) (* @ ) [t ~ y 2n+1 + • • • • (21) 
(2n + ly. 

• The "emergent universe" of [23] can be thought of as an extremality event that has 
been pushed back into the infinite past. 

These definitions have been chosen to match with and simplify the definitions in 
articles [18] and [19]. Note that for bounces these definitions imply that there will 
be some open interval such that 

V* G - A,U) U (*.,*. + A); a(t)>0, (22) 

while for turnarounds there will be some open interval such that 

Vt e (t - A, t ) U (t , t + A); a(t)<0. (23) 

Note that unless the bounce or turnaround is of order one we cannot guarantee that 
at the extremality event itself a(t.) > or a(t ) < 0. For inflexion events we can only 
assert the weaker condition of the existence of some open interval such that 

Vt G (t® - A, t@) U (t®, t® + A); a(t) has fixed sign. (24) 
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Summary: With these definitions in place we can put all of the singular cosmological 
milestones (bangs, crunches, rips, sudden singularities), and the non-singular 
cosmological milestones (extremality events) into a single framework based on 
generalized power series — and we see purely on kinematic grounds that major items of 
interest will be the first two indicial exponents (770 and 771). Indeed 770 > for bangs and 
crunches, r? < for rips, and r? = for sudden singularities (and also, with additional 
conditions, for extremality events). For the nonsingular milestones (bounces, inflexions, 
turnarounds) the parameterization can be somewhat simpler, needing only one positive 
integer — the order of the extremality event, to specify the qualitative behaviour of 
the Taylor series a(t). In the next two sections we will use these parameters to explore 
the generic properties of the cosmological milestones, from the presence (or absence) of 
curvature singularities, to the violation (or otherwise) of the energy conditions. 



3. Kinematics — spacetime curvature 

To analyze the spacetime curvature, let us first consider the Hubble parameter H = a/a. 
From the definition of a generic cosmological milestone 

a(t) = c |t - t©p + a\t - t Q \ m + • • • (25) 

we have (assuming t > t© for simplicity, otherwise one need merely formally reverse the 
flow of time)^f 

a(t) = coVo(t - teY 1 "' 1 + c lVl (t - i©)^ 1 + ... (26) 
Keeping only the most dominant term, we have for r)o ^ 

a co(t - £©)* t - t© 

That is, for bangs, crunches, and rips the Hubble parameter exhibits a generic 1/ (t — t©) 
blow up. For i] = (corresponding to either a sudden singularity or an extremality 
event) we need to go to the next highest term in the numerator (a term which depends 
on r)i, which is guaranteed to be greater than zero by our definitions) to obtain 

H ~ c ^ t ~ t ^ tn ~ 1 = m C -±(t- t Q r~\ ^ = 0; Vl > 0). (28) 

Co Co 

In particular, this guarantees that power law behaviour is completely generic near the 
cosmological milestone, and while the value of the exponent is typically —1, there is 
an exceptional class of milestones (the sudden singularities and extremality events) for 

In fact for explicit calculations in the vicinity of any cosmological milestone it is always possible 
to choose the direction of time to force t > to and so dispense with the need to take the absolute 
value \t — to\, at least for one-sided calculations. This is not a physical restriction on the cosmological 
milestone, just a mathematical convenience which we shall adopt henceforth without further explicit 
discussion. 
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lim H 

t^tQ 



which the exponent will differ. Note that it is not automatic that the Hubble parameter 
diverge at the cosmological milestone. Indeed 

+00 T] > 0; 

sign(ci)cx) i] = 0; rji € (0, 1); 

ci/co Vo = 0; rji = 1; (29) 

770 = 0; 771 > 1; 

—00 770 < 0. 

In particular the Hubble parameter has a finite limit iff 770 = 0, 771 > 1, corresponding 
to a particular subset of the sudden singularities. We have presented this in some 
detail because the style of argument is generic. Consider for instance the cosmological 
acceleration a: 

a(t) = Co?7o ( ?7o -l)|t-t r»- 2 + Cl ? 7l (77 1 -l)|t-t r' 1 - 2 + C2 77 2 (r7 2 - l)|t-t |^" 2 + . . .(30) 

Then provided 770 7^ and 770 7^ 1 
a 770(770 



(31) 



(t-t Q r~ 2 . (32) 



(t-t & r- 2 . (33) 



a (t- t Q f ' 
If 770 = 0, then provided 771 7^ 1 

a 771(771 - l)ci 

a Co 
If both 770 = and 771 = 1 

a 772(772 - l)c 2 

a Co 
Finally, if 770 = 1, then (since 771 > 1) 

a ^ m(m - i) Cl r - 3- (34) 

a c 

Again the generic message is that the behaviour oi a/a near the milestone will be some 
power law, though the precise exponent of that power law will depend on the interplay 
between the various indicial exponents 77,. Note that there is at least one situation 
in which we have to calculate up to the third exponent 772. We can now consider the 
so-called deceleration parameter 

a a 



a 



(35) 



Following the same sort of analysis, for 770 7^ and 770 7^ 1 we have the "generic" result: 
1 — ^70 



too ^0,1). 



(36) 



For the "exceptional" cases we easily see: 
- If 770 = 0, then provided 771 7^ 1 
(771 - l)c 



m 



(37) 
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If both ^o = and rji = 1 

C2C0 
q — 

Finally, if r] 



V2-2 



772(772-1) (t-t Q ) 
1, then (since 771 > 1) 

- ^-^ ( t -t Q r-\ 

c 



(38) 



(39) 



Again we see the ubiquity of power law behaviour, with a "generic" case and a limited 
number of special cases. Note that the deceleration parameter can be either finite or 
infinite as the cosmological milestone is approached: 



lim q 

t-t Q 



(1 -Vo)/Vo 








sign(ci[l - ?7i])oo 


Vo = 


7717^1; 







770 = 


m = i; 


772 > 2; 


-2c 2 c /cl 


Vo = 


vi = i; 


772 = 2; 


-sign(c 2 )oo 


770 = 


771 = 1; 


772G(1,2) 



(40) 



Note that the deceleration parameter has an infinite limit only for a certain subset of 
the sudden singularities 770 = 0. For bangs and crunches the limit is always finite." 1 " 



Polynomial curvature singularities: To now formally decide whether the cosmological 
milestones we have defined (bangs, crunches, rips, sudden singularities, extremality 
events) are curvature singularities in the sense commonly encountered in general 
relativity we will have to look at the Riemann tensor. Because of the symmetries of 
FRW geometry, the Weyl tensor is automatically zero and so it suffices to consider the 
Ricci tensor (and implicitly the Ricci scalar). But because of spherical symmetry, and 
the perhaps less obvious translational symmetry, the only two non-zero orthonormal 
independent components of the Ricci tensor are 



Rt 



and 



Rao — Rj 



(41) 



Thus to test for all possible polynomial curvature singularities it suffices to test for 
singularities in, for instance, R^ and R?f. Alternatively one could consider R^ and the 
Ricci scalar R, or even R^ and G^. Indeed a brief computation reveals that 

a 



Rii - -3 



Gu = 3 



+ 



(42) 
(43) 



In fact these two particular linear combinations of the orthonormal components of the 
Ricci tensor satisfy the useful properties that: 

• they are linearly independent; 

• Rff is independent of the curvature of space (no k dependence); 

+ Of course this is largely because the definition of the deceleration parameter was carefully chosen to 
eliminate the leading t behaviour whenever possible. 
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• Gn * s independent of a (minimizing the number of derivatives involved); 

• testing these two objects for finiteness is sufficient to completely characterize all 
polynomial curvature singularities in a FRW geometry. 

This is observation already enough to tell us that almost all of the cosmological 
milestones defined above are in fact polynomial curvature singularities. Indeed the 
number of situations in which Rn remains finite is rather limited. The generic result is 

lim % = sign(?7o[l-?7o])oo for r?o 7^ 0; t/ ^ 1, (44) 

t^tQ 



which focuses attention on 7/ 



and 7/o = l- From our previous results we see that 



lim R 

t-t 



tt 



—00 


m > i; 





770 = 1; ili > 3; 


-18ci/co 


770 = 1; 7]i = 3; 


— sign(ci)oo 


77 = 1; 77! G (0,3); 


+00 


1oe(0,l); 





% = 0; t/! > 2; 


-6ci/co 


Vo = 0; 7/1 = 2; 


sign(ci)oo 


7/0 = 0; 771 e (1,2); 





% = 0; t/! = 1; 


-6c 2 /c 


?7o = 0; 7/1 = 1; 


-sign(c 2 )oo 


% = 0; 7/1 = 1; 


— sign(ci)oo 


7/o = 0; 77! e (0,1); 


— 00 


7/0 < 0. 



(45) 



% > 2; 
^2 = 2; 

%e(l,2); 



Therefore is finite provided 

• Vo = 0, 7/1 > 2; 

• ?7o = 0, 771 = 1, 7/2 > 2; 

• 7/ = l and 7/! > 3. 

To analyze Gn recall that the condition that the Hubble parameter H remain finite was 
7/0 = 0, 771 > 1. Thus if k > then Gn will remain finite iff H remains finite. In contrast, 
if k < then is a difference of squares and one term can be balanced against the 
other. Indeed near the milestone 

k 



while 



Gii ~ 



G B ~ 



+ 



(t-t Q y ' ci(t-t Q y™ 

riic\{t - t e ) 2 ^~V + k 



7/o = 0. 



(46) 



(47) 
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After a brief calculation 



sign(A;)oo 


Vo > l; 


fc^O; 




+00 


^0 > l; 


fc = 0; 




signfco + k)oo 


?7o = l; 


9 1 7 / r\ 

c 2 + k^0; 







Vo = l; 


eg + fc = 0; 


771 > 3; 


18ci 


?7o = l; 


eg + fc = 0; 


771 = 3; 


sign(ci)oo 


?7o = l; 


eg + k = 0; 


Vi e (1,3); 


+00 


%e(0,l); 






3k/c 2 


^o = 0; 


771 > 1; 




3(c? + k)/cl 


^o = 0; 


77i = l. 




+00 


^o = 0; 


771 e (0,1). 





(48) 



Thus the necessary and sufficient conditions for to remain finite at the cosmological 
milestone are: 

• Vo = 0,Vi> i; 

• i] = 1, eg + fc = 0, and 77! > 3. 

(Note that eg + k = implies k = —1, and c = 1.) 

Combining with the previous result, the necessary and sufficient conditions for both G% 
and Rj$ to remain finite, so that a cosmological milestone is not a polynomial curvature 
singularity, are: 

• Vo = 0, 771 > 2; 

• Vo = 0, 7/! = 1, 772 > 2; 

• 770 = 1, k = —1, Co = 1, and 771 > 3. 

This can only happen for a particular sub-class of sudden singularities (770 = 0), and for a 
rather exceptional type of big bang/crunch, one that asymptotes to a Milne universe (see 
discussion below). However, once one generalizes the concept of polynomial curvature 
singularity to that of a derivative curvature singularity, most of these exceptional cases 
will be excluded. 



Derivative curvature singularities: A derivative curvature singularity is defined by some 
polynomial constructed from finite-order derivatives of the curvature tensor blowing up. 
Because of the symmetries of the FRW universe the only interesting derivatives will be 
time derivatives, and so the only objects we need consider are 

and (49) 



But 



while 



d n t d n t 



d n R- a < - n+2 ' ) 

= —3 h (lower-order derivatives); (50) 



d n G u Jo( n+1 » 



d n t a 



2 



+ (lower-order derivatives). (51) 
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Thus to avoid a nth-order derivative curvature singularity we must at the very least keep 
a("+ 2 )/a finite, and furthermore all related lower-order derivatives of the form a^/a, 
with j < n + 2, must also be finite. To prevent any arbitrary-order derivative singularity 
from occurring, we must force all /a to remain finite. This condition holds in addition 
to the constraint coming from polynomial curvature singularities derived above. 

Since the constraint from polynomial curvature singularities had already forced 
770 = or r]o — 1, then this implies that we must force all to remain finite, and so 
must force all the indicial exponents rji to be non-negative integers, thus making a(t) a 
Taylor series. This condition holds in addition to the constraint coming from polynomial 
curvature singularities derived above, with the consequence that almost all cosmological 
milestones are indeed curvature singularities. Indeed, inserting the condition that a(t) 
be representable by a Taylor series into and shows that the only two situations 
in which a cosmological milestone is not a derivative curvature singularity is if: 

• ?7o = 0, rji G Z + ; corresponding to an extremality event (bounce, turnaround, or 
inflexion event) rather than a bang, crunch, rip, or sudden singularity; 

• rjo — 1, k — — 1, Cq = 1, rji E Z + , and rji > 3; corresponding to a FRW geometry 
that smoothly asymptotes near the cosmological milestone to the Riemann-flat 
Milne universe. (The special case a(t) — t, k — — 1 is called the Milne universe and 
is actually a disguised portion of Minkowski space. It corresponds to the interior of 
the future light cone based at some randomly specified point in Minkowski space, 
with a spatial foliation defined by the proper time hyperboloids based on that event. 
Since this spacetime is flat, it corresponds to a universe which on the largest scales 
is empty. This is not a popular cosmological model.) 

In other words, almost all cosmological milestones are physical singularities, apart from a 
very limited sub-class corresponding to either extremality events (bounces, turnarounds, 
inflexion events) or an asymptotically empty universe. It is important to realise that up 
to this stage our considerations have been purely kinematic — up to this point we have 
not used or needed the Friedmann equations or the Einstein equations. 

4. Dynamics — energy conditions 

To now start to include dynamics we relate the geometry to the density and pressure 
using the Friedmann equations and then ask what happens to the energy conditions at 
the cosmological milestones. Even though we may have good reason to suspect that 
the energy conditions are not truly fundamental [24] - they make a very good first 
pass at the problem of quantifying just how "strange" physics gets at the cosmological 
milestone. The standard energy conditions are the null, weak, strong, and dominant 
energy conditions which for a FRW spacetime specialise to [18, 19, 24]: 

• [NEC] p + p > 0. 

In view of the Friedmann equations this reduces to 
a? + k 

a < ; that is k > a a — a 2 . (52) 

a 
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[WEC] This specialises to the NEC plus p > 0. 
This reduces to the NEC plus the condition 

a 2 + k > 0; that is k > —a 2 . 

[SEC] This specialises to the NEC plus p + 3p > 0. 
This reduces to the NEC plus the deceleration condition 

a < 0. 

[DEC] p±p > 0. 

This reduces to the NEC plus the condition 

2{a 2 + k) 



(53) 



(54) 



2(a 2 + k) , . , (aa + 2a 2 ) . . 

a > i L that is k > -- -. (55) 

a 2 

Of these, it is the NEC that is the most interesting (because this is the weakest of the 

standard energy conditions it leads to the strongest theorems), and we shall analyze 

this condition in most detail. 



NEC: Now near any generic cosmological milestone 

while in the degenerate cases 

a a - a 2 ~ c c 1 7/1(771 - 1) (f - t Q ) {m ~ 2) ; (t/ = 0; 771 ^ 1), 



and 



Therefore 



a a - a 2 ~ c c 2 772(772 - 1) (t - t Q )^~ 2) - cf; 



(56) 
(57) 

(7fo = 0; 771 = 1). (58) 



lim (a a — a 2 ) 






Vo > 


i; 








-cl 


Vo = 


i; 








—00 


Vo e 


(0,1); 











^0 = 


0; 


Vi > 


2; 




2c ci 


^0 = 


0; 


Vi = 


2; 




+sign(ci)oo 


^0 = 


0; 


771 G 


(1,2) 




-cl 


^0 = 


0; 


»7i = 


1; 


?72>2; 


2CqC2 — c x 


Vo = 


0; 


»7i = 


1; 


r? 2 = 2; 


+sign(c 2 )oo 


Vo = 


0; 


Vi = 


1; 


r7 2 e(l,2 


— sign(ci)oo 


Vo = 


0; 


Vi e 


(0,1) 




+00 


Vo < 


0. 









(59) 



Hence the NEC is definitely satisfied (meaning the inequality is strict) at a generic 
cosmological milestone iff: 

• 770 > 1, k = +1; 

• 7/0=1: for = +1 and k — 0, and also for fc = — 1 with the proviso that c > 1; 

• Vo £ (0, 1): for any value of fc; 
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• r] = subject to the additional constraints: 

— 771 > 2 and k = +1; 

— rji — 2 and > 2coCi; 

— 771 G (1, 2) and ci < 0; 

— 771 = 1, 772 > 2: for /c = +1, A; = 0, and for k = —1 with the proviso |ci| > 1; 

— 771 = 1, 772 = 2: for fc > 2c c 2 - cf; 

— 771 = 1, 772 e (1, 2), and c 2 < 0; 

— i]i G (0, 1), and Ci > 0. 

Hence the NEC is marginally satisfied (meaning the non-strict inequality is actually an 
equality) at a generic cosmological milestone iff: 

• 770 > 1, k — 0; 

• 770 = 1, k — — 1 with the proviso that c = 1; 

• ?7o = subject to the additional constraints: 

— rji > 2 and /c = 0; 

— ?7i = 2 and k = 2cqC\ (which requires k 7^ 0); 

— 771 = 1, 7/2 > 2: for k = —1 and ci = ±1; 
-7/i = l,7/ 2 = 2, and fc = 2c c 2 - cf. 

In all other situations the NEC is definitely violated. To make this explicit, the NEC is 
definitely violated (the inequality is strictly violated) at a generic milestone iff: 

• 770 > 1, k = -1; 

• 7/o = l: for k = — 1 with the proviso that c < 1; 

• 7/o = subject to the additional constraints: 

— 7/1 > 2 and k = — 1; 
-7/i = 2 and A; < 2c Ci; 

— 7/1 G (1, 2) and c 1 < 0; 

— 7/1 = 1, 7/2 > 2: for k = — 1 with the proviso |ci| < 1; 

— 771 = 1, 7/2 = 2: for fc < 2c c 2 - cj; 

— 771 = 1, 772 G (1, 2), and c 2 > 0; 

— 7/1 G (0, 1), and c x < 0. 

• 7/0 < 0. 

This pattern is in agreement with the standard folklore, and is a systematic expression 
of results that are otherwise scattered throughout the literature. For instance, it is 
immediate that "big rips" always violate the NEC (and hence all the other energy 
conditions), and that for big bangs and crunches the range 770 G (0, 1) is preferred, since 
it is only in this range that the NEC holds independent of the spatial curvature. Note 
that for sufficiently violent big bangs (7/0 > 1) and hyperbolic spatial curvature the 
NEC is violated — this indicates that "phantom matter" need not always lead to a "big 
rip" [7, 8]; it might in contrast lead to a particularly violent bang or crunch. It is for 
7/o = (which corresponds to either sudden singularities or extremality events) that the 
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analysis becomes somewhat tedious. Certainly, many (though not all) of the sudden 
singularities violate the NEC. 



WEC: If in addition we consider the WEC, we need to look at 

while in the degenerate case 

-d 2 --^-^" 1 ); (770 = 0). 

Therefore 



lim {—a ) 






Vo > l; 




-r 2 


= l; 




— oo 


%e (0,1); 







Vo = 0; 


m > l; 


-r 2 


Vo = 0; 


Vi = !; 


— oo 


Vo = 0; 


vi e (o, l) 


— oo 


Vo < 0. 





(60) 
(61) 



(62) 



Since we want the condition k > —a 2 to hold in addition to the NEC, it follows that the 
WEC is equivalent to the NEC for bangs crunches and rips, and that the only changes 
arise when dealing with 770 — (sudden singularities or extremality events). Certainly, 
many (though not all) of the sudden singularities violate the WEC. 



SEC: When analyzing the SEC, it is sufficient to note that near any generic 
cosmological milestone: 

d~ 770(770-1) c 1<"°- 2 >; (770^0, 770^1), (63) 

while in the degenerate cases 



a ~ 77 1 (77 1 - 1) d (t-t G )^- 2) ; 
a ~ 772(772 - 1) c 2 (t-t )te- 2 ) ; 



Vo = 0; 7/x ^ 1; 
Vo = 1; 

(?7o = 0; 771 = 1;), 



Therefore 



lim [sign(d)] = < 



( +1 

sign(ci) 
-1 

sign (771 [771 - l]ci) 
sign(c 2 ) 
+ 1 



(64) 
(65) 



(66) 



Vo > 1 

Vo = l; 

Voe(0,l): 
Vo = 0; 771 ^ 1; 
?7o = 0; 771 = 1; 
770 < 0. 

Since we want the condition a < to hold in addition to the NEC, it follows that the 
SEC is definitely violated for all "rips" , and is also definitely violated for all "violent" 
bangs and crunches with 770 > 1. The SEC is definitely satisfied for bangs and crunches 
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in the range r] G (0, 1). For turnarounds or bounces rji = 2n, n G Z + , and the sign of 
Ci governs possible SEC violations: Turnarounds satisfy the SEC while bounces violate 
the SEC. (See also [18, 19].) For inflexion events, where rji — 2n + 1, n G Z + , the SEC 
is violated either just before or just after the inflexion event. (See also [18, 19].) 



DEC: Finally, when analyzing the DEC, it is sufficient to note that near any generic 
cosmological milestone: 



ad + 2d 2 ) Voi^Vo - 1) c o 



2 2 
while in the various degenerate cases 

(ad + 2d 2 ) 771(771 - 1) c c x 



{t _ t&) 2(r ]0 -l). (??o ^ 0) ^1/3), 



(t - t Q) 



{ad + 2d 2 ) 772(772 ~ 1) c c 2 ( 2) 2 
2 ~ ^ (* ~ *©) 



(t/o = 0; 771 ^ 1), 

(770 = 0; 771 = 1), 



(ad + 2d 2 ) (37/ 1 + 2)(377 1 -l)coc 1 



Therefore 
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(t-t ) (m - 5/3) ; (770 = 1/3). 



lim 



(ad + 2d 2 






% > 1 




-c 2 


77o = 1 




—00 


770 e (1/3,1); 





770 = 1/3; 77! > 5/3 


14 

--9-C0C1 


770 = 1/3; 771 = 5/3 


— sign (ci) 00 


770 = 1/3; 77x < 5/3 


+00 


77o e (0,1/3); 





77o = 


771 > 2; 


-C Ci 


77o = 


771 =2; 


— sign (ci) 00 


77o = 


771 e (1,2) 




77o = 


771 = 1; 


-c c 2 - c 2 


77o = 


771 = i; 


-sign(c 2 )oo 


77o = 


771 = 1; 


+sign(ci)oo 


77o = 


771 e (0,1); 


—00 


770 < 





(67) 

(68) 
(69) 
(70) 



(71) 



772 > 2; 
772 = 2; 
772 e (1,2); 



Remember that to satisfy the DEC one needs to satisfy the NEC in addition to the 
constraint coming from the above. Let us define 

(ad + 2d 2 ) 



K = max < lim lad — d 2 ] , lim 

1 t^t® 1 J t^t© 



(72) 



Then satisfying the DEC at the cosmological milestone is equivalent to the constraint 
k > K. (73) 
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For bangs and crunches we calculate: 



K= { 



— oo 


14 „ _ 
— -$CqCi 

— sign(ci)oo 
+00 



770 > 1; 

-co % = l; 

770 g (1/3,1); 

r]o = 1/3; 77! > 5/3; (74) 
770 = 1/3; 7/! = 5/3; 
770 = 1/3; 771 < 5/3; 

770 g (0,1/3). 

In particular, if one wishes to use "normal" matter (that is, matter satisfying all of the 
energy conditions) to drive a bang or a crunch independent of the sign of space curvature 
k, then one is forced in a model independent manner into the range 770 G (1/3, 1). 

— For rips we have: 

K = +oo 770 < 0. (75) 

Note that rips, because they violate the NEC, always violate the DEC, so the particular 
constraint derived above is not the controlling feature. 

— Finally, for sudden singularities and extremality events: 



K 



Co max{2ci, — c±} 



-00 



-1 



c max{2ci, — c±} 

+00 

+00 



V2 > 2; 

m e (1,2); 



770 = 0; 77! > 2; 

Vo = 0; 771 = 2; 

770 = 0; (1,2) 

770 = 0; 77! = 1; r] 2 > 2; (76) 

1 ^o = 0; 771 = 1; 

Vo = 0; 7/! = 1; 

770 = 0; 771 G (0,1). 

As we have now seen (several times), the 770 = case is "special" and requires extra care 
and delicacy in the analysis — this is the underlying reason why "sudden singularities" 
are so "fragile" , and so dependent on the specific details of the particular model. Indeed 
there are several classes of cosmological milestone for which the DEC is satisfied. For 
instance, a complete catalogue of the bangs and crunches for which the DEC is satisfied 
is: 

• 770 > 1, k = 0, +1. 

• 770 = 1, k — 0, +1, and k = — 1 if Co > 1. 

• 770 G (1/3,1), for all k. 

• 770 = 1/3, 77! > 5/3, k = 0,+l. 

• 770 = 1/3, 771 = 5/3, k > -^c ci. 

• 770 = 1/3, 7/! G (1/3,5/3), ci > 0, for all k. 

In particular, if the DEC is to hold independent of the sign of space curvature k all the 
way down to the singularity, then in the vicinity of the singularity the dominant term 
in the scale factor a(t) is bounded from both above and below by 

Co(t ~ t*) 1/3 < OdominantW < C (t - Q. (77) 
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Similarly, a complete catalogue of the sudden singularities for which the DEC is satisfied 
is: 

• 770 = 0, 771 > 2, k = 0,+l. 

• rjo = 0, rji = 2, with k > Co max{2ci, — ci}. 

(In particular this requires = +1 though that is not sufficient.) 

• 770 = 0, 771 = 1, 772 > 2, k — 0, +1, and k — — 1 if |ci| > 1. 

• r/o = 0, r/i = 1, i] 2 = 2, fc > c max{2ci, -ci} - c\. 

This appears at first glance, to contradict the analysis of Lake [15], who claimed that all 
sudden singularities violate the DEC. The resolution of this apparent contradiction lies 
in the question "just how sudden is the sudden singularity?" Lake takes as his definition 
a(t — > to) = —00, corresponding in our analysis to r] = with either 771 G (0, 1) U (1, 2) 
or rji = 1 with 772 < 2. In this situation our results certainly agree with those of Lake - 
if d(t — > t & ) = —00 then the DEC is certainly violated. However, our "counterexamples" 
where the DEC is satisfied all satisfy a(t — > t Q ) = finite, and for these counterexamples it 
is only some higher derivative a^ a \t — > t & ) for n > 3 that diverges. Sudden singularities 
of this type are sufficiently "gentle" that at least for some of them even the DEC can 
be satisfied all the way to the singularity - - a conclusion completely in agreement 
with Barrow and Tsagas [13], though now we have a complete characterization of those 
situations for which the DEC can be satisfied.* 

5. Discussion 

In this article we have explored three issues: First we have developed a complete 
catalogue of the various "cosmological milestones" in terms of a generalized power 
series expansion of the FRW scale factor in the vicinity of the milestone. This power 
series expansion is sufficiently general to accommodate all commonly occurring models 
considered in the literature. 

Second, with the notion of a generalized power series in hand, it is possible at 
a purely kinematic level to address the question of when a "cosmological milestone" 
corresponds to a curvature singularity, and what type of singularity is implied. In this 
way is is possible to classify all "cosmlogical milestones" as to whether they are singular 
or not, with the class non-singular milestones being strictly limited (to extremality 
events or asymptotically Milne bangs/crunches). 

Third, this definition of cosmological milestones in terms of generalized power series 
enables us to perform a complete model-independent check on the validity or otherwise 
of the classical energy conditions. In particular we provide a complete catalogue of 
those bangs/crunches, sudden singularities and extremality events for which the DEC 
is satisfied, and a complete catalogue of those bangs/crunches, sudden singularities and 
extremality events for which the NEC is satisfied. Depending on one's attitude towards 

* See also the discussion by Sahni and Shtanov [10] in a braneworld context. 
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the energy conditions [24], one could use this catalogue as a guide towards deciding on 
potentially interesting scenarios to investigate. 

Finally, we remind the reader of what is perhaps the most general lesson to be 
learned: If in the vicinity of any cosmological milestone, the input scale factor a(t) is a 
generalized power series, then all physical observables (H, q, the Riemann tensor, etc.) 
will likewise be generalized power series, with related indicial exponents that can be 
calculated from the indicial exponents of the scale factor. Whether or not the particular 
physical observable then diverges at the cosmological milestone is "simply" a matter of 
calculating its dominant indicial exponent in terms of those occuring in the scale factor. 

Acknowledgements 

This research was supported by the Marsden Fund administered by the Royal Society 
of New Zealand. 

References 

[1] C. W. Misner, K. S. Thome, and J. A. Wheeler, Gravitation, (Freeman, San Francisco, 1972), 
[2] R. M. Wald, General relativity, (Chicago University Press, 1984). 

[3] S. M. Carroll, Spacetime and geometry: An introduction to general relativity, (Addison- Wesley, 
San Francisco, 2004). 

[4] J. B. Hartle, Gravity: An introduction to Einstein's general relativity, (Addison- Wesley, San 
Francisco, 2003). 

[5] M. Visser, "Jerk and the cosmological equation of state," Class. Quant. Grav. 21 (2004) 2603 
[arXiv:gr-qc/0309109]. 

[6] M. Visser, "Cosmography: Cosmology without the Einstein equations," arXiv:gr-qc/0411131. 
[7] R. R. Caldwell, "A Phantom Menace?," Phys. Lett. B 545 (2002) 23 [arXiv:astro-ph/9908168]. 
[8] R. R. Caldwell, M. Kamionkowski and N. N. Weinberg, "Phantom Energy and Cosmic Doomsday," 

Phys. Rev. Lett. 91 (2003) 071301 [arXiv:astro-ph/0302506]. 
[9] P. F. Gonzalez-Diaz, "You need not be afraid of phantom energy," Phys. Rev. D 68 (2003) 021303 

[arXiv: astro-ph/0305559] . 

M. Bouhmadi-Lopez and J. A. Jimenez Madrid, "Escaping the Big Rip?," JCAP 0505 (2005) 
005 [arXiv:astro-ph/0404540]. 

L. P. Chimento and R. Lazkoz, "On big rip singularities," Mod. Phys. Lett. A 19 (2004) 2479 
[arXiv:gr-qc/0405020]. 

P. H. Frampton and T. Takahashi, "Bigger rip with no dark energy," Astropart. Phys. 22 (2004) 
307 [arXiv:astro-ph/0405333]. 

S. Nojiri and S. D. Odintsov, "Quantum escape of sudden future singularity," Phys. Lett. B 595 
(2004) 1 [arXiv:hep-th/0405078]. 

P. H. Frampton, "Enhanced big rip without conventional dark energy," arXiv:astro-ph/0407353. 
E. Babichev, V. Dokuchaev and Y. Eroshenko, "Dark energy cosmology with generalized linear 
equation of state," Class. Quant. Grav. 22 (2005) 143 [arXiv:astro-ph/0407190]. 
H. Caldcron and W. A. Hiscock, "Quantum fields and 'Big Rip' expansion singularities," Class. 
Quant. Grav. 22 (2005) L23 [arXiv:gr-qc/0411134]. 

S. Nojiri, S. D. Odintsov and S. Tsujikawa, "Properties of singularities in (phantom) dark energy 
universe," Phys. Rev. D 71 (2005) 063004 [arXiv:hcp-th/0501025]. 

R. Curbelo, T. Gonzalez and I. Quiros, "Interacting phantom energy and avoidance of the big 
rip singularity," arXiv:astro-ph/0502141. 



Necessary and sufficient conditions for big bangs, bounces, crunches, rips... 19 

S. Nojiri and S. D. Odintsov, "Inhomogeneous equation of state of the universe: Phantom era, 
future singularity and crossing the phantom barrier," Phys. Rev. D 72 (2005) 023003 [arXiv:hep- 
th/0505215]. 

[10] Y. Shtanov and V. Sahni, "Unusual cosmological singularities in braneworld models," Class. Quant. 

Grav. 19 (2002) L101 [arXiv:gr-qc/0204040]. 
[11] J. D. Barrow, "Sudden Future Singularities," Class. Quant. Grav. 21 (2004) L79 [arXiv:gr- 

qc/0403084]. 

[12] J. D. Barrow, "More general sudden singularities," Class. Quant. Grav. 21 (2004) 5619 [arXiv:gr- 
qc/0409062]. 

[13] J. D. Barrow and C. G. Tsagas, "New Isotropic and Anisotropic Sudden Singularities," Class. 

Quant. Grav. 22 (2005) 1563 [arXiv:gr-qc/0411045]. 
[14] S. Nojiri and S. D. Odintsov, "Quantum escape of sudden future singularity," Phys. Lett. B 595 

(2004) 1 [arXiv:hep-th/0405078]. 

L. Fernandez-Jambrina and R. Lazkoz, "Geodesic behaviour of sudden future singularities," 
Phys. Rev. D 70 (2004) 121503 [arXiv:gr-qc/0410124]. 

H. Stefancic, "'Expansion' around the vacuum equation of state: Sudden future singularities 
and asymptotic behavior," Phys. Rev. D 71 (2005) 084024 [arXiv:astro-ph/0411630]. 
[15] K. Lake, "Sudden future singularities in FLRW cosmologies," Class. Quant. Grav. 21 (2004) L129 
[arXiv:gr-qc/0407107]. 

[16] M. Visser and N. Yunes, "Power laws, scale invariance, and generalized Frobenius series: 
Applications to Newtonian and TOV stars near criticality," Int. J. Mod. Phys. A 18 (2003) 
3433 [arXiv:gr-qc/0211001]. 

[17] J. D. Bekenstein, "Nonsingular General Relativistic Cosmologies," Phys. Rev. D 11 (1975) 2072. 

[18] C. Molina-Paris and M. Visser, "Minimal conditions for the creation of a Friedman-Robertson- 
Walker universe from a 'bounce'," Phys. Lett. B 455 (1999) 90 [arXiv:gr-qc/9810023]. 

[19] D. Hochberg, C. Molina-Paris and M. Visser, "Tolman wormholes violate the strong energy 
condition," Phys. Rev. D 59 (1999) 044011 [arXiv:gr-qc/9810029]. 

[20] L. Parker and Y. Wang, "Avoidance Of Singularities In Relativity Through Two-Body 
Interactions," Phys. Rev. D 42, 1877 (1990). 

R. Brustein, S. Foffa and A. E. Mayo, "Causal entropy bound for non-singular cosmologies," 
Phys. Rev. D 65 (2002) 024004 [arXiv:hep-th/0108098]. 

J. C. Hwang and H. Noh, "Non-singular big-bounces and evolution of linear fluctuations," Phys. 
Rev. D 65 (2002) 124010 [arXiv:astro-ph/0112079]. 

C. Gordon and N. Turok, "Cosmological perturbations through a general relativistic bounce," 
Phys. Rev. D 67 (2003) 123508 [arXiv:hcp-th/0206138]. 

S. Tsujikawa, R. Brandenberger and F. Finelli, "On the construction of nonsingular pre-big-bang 
and ekpyrotic cosmologies and the resulting density perturbations," Phys. Rev. D 66, 083513 
(2002) [arXiv:hcp-th/0207228]. 

Y. Shtanov and V. Sahni, "Bouncing braneworlds," Phys. Lett. B 557 (2003) 1 [arXiv:gr- 
qc/0208047]. 

A. J. M. Medved, "Anatomy of a bounce," Class. Quant. Grav. 21 (2004) 2749 [arXiv:hcp- 
th/0307258]. 

G. Vcncziano, "A model for the big bounce," JCAP 0403 (2004) 004 [arXiv:hcp-th/0312182]. 
M. Gasperini, M. Giovannini and G. Veneziano, "Cosmological perturbations across a curvature 
bounce," Nucl. Phys. B 694 (2004) 206 [arXiv:hep-th/0401112]. 

L. E. Allen and D. Wands, "Cosmological perturbations through a simple bounce," Phys. Rev. 
D 70 (2004) 063515 [arXiv:astro-ph/0404441]. 

M. G. Brown, K. Freese and W. H. Kinney, "The phantom bounce: A new oscillating cosmology," 
arXiv:astro-ph/0405353. 

B. K. Tippett and K. Lake, "Energy conditions and a bounce in FLRW cosmologies," arXiv:gr- 
qc/0409088. 



Necessary and sufficient conditions for big bangs, bounces, crunches, rips. . . 20 

L. E. Allen, "Cosmo-logical perturbations through a simple bounce," AIP Conf. Proc. 736 (2005) 
182. 

[21] L. H. Ford, "Unstable Fields And The Recollapse Of An Open Universe," Phys. Lett. A 110 (1985) 
21. 

J. D. Barrow, "The Premature Recollapse Problem In Closed Inflationary Universes," Nucl. 
Phys. B 296 (1988) 697. 

G. A. Burnett, "Recollapse of the closed Tolman space-times," Phys. Rev. D 48 (1993) 5688 
[arXiv:gr-qc/9308003]. 

J. Miritzis, "The recollapse problem of closed FRW models in higher-order gravity theories," 
arXiv:gr-qc/0505139. 

[22] V. Sahni, H. Fcldman and A. Stebbins, "Loitering universe," Astrophys. J. 385 (1992) 1. 

V. Sahni and Y. Shtanov, "Did the Universe loiter at high redshifts ?," Phys. Rev. D 71 (2005) 

084018 [arXiv:astro-ph/0410221]. 
[23] G. F. R. Ellis and R. Maartens, "The emergent universe: Inflationary cosmology with no 

singularity," Class. Quant. Grav. 21 (2004) 223 [arXiv:gr-qc/0211082]. 

G. F. R. Ellis, J. Murugan and C. G. Tsagas, "The emergent universe: An explicit construction," 
Class. Quant. Grav. 21 (2004) 233 [arXiv:gr-qc/0307112]. 
[24] C. Barcclo and M. Visser, "Twilight for the energy conditions?," Int. J. Mod. Phys. D 11 (2002) 
1553 [arXiv:gr-qc/0205066]. 



